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Abstract 

We extend the classical risk minimization model with scalar risk measures to the general 
case of set-valued risk measures. The problem we obtain is a set-valued optimization model 
and we propose a goal programming-based approach with satisfaction function to obtain a 
solution which represents the best compromise between goals and the achievement levels. 
Numerical examples are provided to illustrate how the method works in practical situations. 

Keywords: Risk measure, multi-criteria portfolio optimization, goal programming, satisfaction 
function. 

1 Introduction 

Risk measures are real valued functionals denned on a space of random variables which encloses 
every possible financial position. It may seem naive to use a single number to describe the com- 
plexity of the distributions characterizing those random variables. On the other hand this appears 
as the only way to succeed in the assessment of the capital requirement needed to a bank to 
recover a high possible loss due to risky investments. Without any doubt it is a critical point the 
choice of the axioms defining the risk measures; these have been vividly discussed since the very 
beginning of this theory, opening a broad new branch of research which still triggers the interest 
of the Mathematical Finance world. 

Seminal contributions to this topics were surprisingly given by Bernoulli in 1738 [13| . who per- 
ceived the role of the risk aversion in decision making. This was the starting point of the notion 
of Expected Utility. Later on, in a financial environment, many risk procedure were introduced. It 
was the case of the Mean-Variance criterion (Markovitz, 1952, [35]), the Sharpe's ratio (1964, [35]) 
and the Value at Risk (V@R), defined through the quantiles of a given distribution with a prede- 
fined level of probability. This last method is the most employed in credit institutes and has been 
pointed out as the reference parameter by the Basel Committee on Banking Supervision (Basel II 
2006). 

At the end of the Nineties, Artzner, Delbaen, Eber and Heath produced a rigorous axiomatic 
formalization of coherent risk measures, led by normative intent. The regulating agencies asked 
for computational methods to estimate the capital requirements, exceeding the unmistakable lacks 
showed by the extremely popular V@R. Given a vector space of random variables L, the definition 
of a coherent risk measure p : L — > K requires four main hypotheses to be satisfied: monotonicity, 
cash additivity, positive homogeneity sublinearity (see Section |2j. The relevant role of the axioms 
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was deeply discussed in many papers: Follmer and Schied (2002, [17]), Frittelli and Rosazza Gianin 
(2002, [T5]) independently studied the convex case weakening positive homogeneity and sublinear- 
ity El Karoui and Ravanelli relaxed the cash additivity axiom to cash subadditivity (2009, 15 ) 
when the market presents illiquidity; Maccheroni et al. (2010, [2]) showed how quasiconvexity 
describes better than convexity the principle of diversification, whenever cash additivity does not 
hold. 

Finally two important recent generalizations were introduced by Jouini et al. (2004,(55]), who 
defined set- valued coherent risk measures, and by Hamel and Heyde (2010,(12]) who introduced 
the notion of set-valued convex risk measure. This approach is absolutely natural as far as the 
risk is expressed and hedged in different currencies. 

Diversification plays a crucial role in insurance and financial business and the interpretation of 
this notion was the source of this vivid debate. An agent who considers a fixed basket of financial 
instruments X — (X\, Xd) € tries to reallocate his wealth, by means of a diversified strategy 
a G K d , in order to minimize the risk of his portfolio. Namely, given a real valued risk measure 
(as introduced in [10]) p : L°° ->lwe have the following optimization model 

min Ip (a ■ X) | a 6lf : = lj (1) 

where a • X is the usual scalar product in E d . 

The optimal risk allocation is a classical problem in mathematical economics and it is interesting 
from both practical and theoretical perspectives. In more recent years this problem has also 
been studied in many other contexts such as risk exchange, assignment of liabilities to daughter 
companies, individual hedging problems (see, for instance, the papers by Heath and Ku (2004, [21]), 
Barrieu and El Karoui (2005,[Tl]), Burgert and Ruschendorf (2006,[12 ), Jouini et al. (2007, [23]). 
Acciaio (2007,51)). 

The aim of this paper is to provide a computational procedure, based on the Goal Programming 
(GP) model, to problem ([!]) if the agent has to find the best compromise among different beliefs 
which may come either from the uncertainty on the probabilistic model P, or from different opinion 
that the agent has to face in his institution. These multi-criteria will be aggregated in a unique 
measure of risk which will be described by a set- valued map R : If£ =4 lR n . In particular d stands 
for the number of financial instruments considered, whereas n is the number of different criteria 
(which in general might be larger that d as shown in Section [5]). Thus the interpretation we are 
giving to R appears pretty different to the original one provided in |22j . 

The paper is organized as follows: in Section [2] we recall some basic notions on set valued risk 
measures, Section [3] is devoted to the extension of the optimization problem (JlJ to the case of 
set- valued risk measures, Section [4] introduces a GP model with satisfaction function and finally 
Section [5] presents some numerical results. 

2 Risk Measures 

Let (f2, J 7 , P) a probability space and we denote with L° =: L°(fi, T, P) the space of T measurable 
random variables that are P almost surely finite. We denote by L°° =: L°°(fi, T, P) the space of 
P-almost surely bounded random variables which becomes a Banach lattice once endowed with the 
P-almost sure pointwise partial order and the usual norm of the supremum. In this probabilistic 
framework we recall the definition of risk measure. 

Definition 1 A risk measure is a functional p : L°° —¥ M which satisfies 

i) monotonicity, i.e. X\ < X2 implies p(X\) > p{X2) for every X\,X% G L°° , 
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Moreover a risk measure may satisfy 

ii) convexity, i.e. p(tXi + (1 - t)X 2 ) < tp{X x ) + (1 - t)p(X 2 ) for all t £ [0, 1]. 
in) cash additivity, i.e. p(X + c) = p(X) — c, 
iv) positive homogeneity, i.e. for every a > 0, p(ptX) = ap(X), 
v) sublinearity, i.e. p(X + Y) < p(X) + p(Y). 

Monotonicity represents the minimal requirement for a risk measure to model the preferences 
of a rational agent. If in addition conditions iii), iv) and v) hold then the risk measure is called 
coherent and it is automatically convex. Unfortunately both axioms iv) and v) appear to be 
restrictive and unrealistic: the former does not sense the presence of liquidity risks, the latter does 
not describe the real intuition hidden behind the diversification process. For this reason in most 
of the literature iv) and v) are substituted by ii) which has a natural interpretation: the risk of 
the diversified aggregated position tX\ + (1 — t)X 2 is surely smaller than the combination of the 
two single risks. Cash additivity is a key property which allows to characterize the risk procedure 
in terms of capital requirements 

p(X) = inf{a \X + aeA} 

where A = {X £ L°° \ p(X) < 0} is the acceptance sets. The risk of a position X is thus the 
minimal amount of money that I have to save today in order to make the position acceptable with 
respect to a precise criterion (represented by the acceptance set A) which is usually imposed by 
regulation agencies (External Risk Measures). On the other hand an institution may have some 
specific criteria that need to be enclosed in their model as in the case of Internal Risk Measures. 

In literature, some extensions of the notion of risk measure have been considered in order to 
better describe the complexity of the risk process. We now recall the notion of set valued risk 
measures presented in [22] . Given any subset A C R d we shall denote by L P AA) the collection of 
A- valued random variable X = (X\, Xd) with finite LP norm (or equivalently X; £ LP for every 
i = l, ...,d). Whenever no confusion arises we denote L p d =: L d (R). Notice that for p = +oo we 
end up with essentially bounded random vectors of dimension d. 

In this paper we take into account the theoretical framework developed in [22] and [19] : consider 
a closed convex cone K d C R d (resp. K n C R n ) such that R^ C K d (resp. W% C K n ) and define 
the partial ordering =4d on R d by x =4d iff x € K d (similarly for =^„ on R"). This ordering can 
be naturally extended to in the following way: 

X 4 Y X -Y £ K d P-almost surely 

Hence L^°(X d ) is a cone that consists in all the non-negative random variables in the sense of 
Moreover for any A, B C R n we may define the partial order =<;„ as 

A 4n B & B CA + K n 

We will indicate by (1,1,. ..,1) =: Id E R d (l n G R n ) the vector such that each component is 
equal to 1, by l 3 d € R d (1^ G R d ) the vector such that the j*' 1 component is 1 and the other are 0. 
Finally the sum among sets is to be intended as the usual Minkwoski sum. 

Definition 2 A (d, n)-risk measure is a set valued map R : =^ 1" satisfying the following 
axioms: 

i) for all X e , R{X) is closed and G R(Q) ^ R n : 

ii) for all X,Y £ Lf: X ^ d Y V-almost surely implies R(Y) =<!„ R(X). 
In particular a {d,n)-risk measure is convex if 
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Hi) for all Xi,X 3 £ Lf, Wt G [0, 1], R(tXi + (1 - t)X 2 ) =*„ tiJ(JTi) + (1 - t)R{X 2 ) 
A (d,n)-risk measure is cash additive if 

iv) for all x £ R, j = 1, d and X £ Lf we have R(X + xl j d ) = -xl n + R{X); 
and coherent if 

v) it is sublinear: for all X, Y £ Lf, R(X + Y) 4 n R(X) + R(Y) and 

vi) positively homogenous: for all t > and X G Lf , R(tX) = tR(X). 
Remark 3 We now illustrate the financial meaning of the previous axioms. 

1. Monotonicity implies that for X =^<j Y the efficient points of R(X) 'dominate 'those of R(Y). 
This domination might be described as follows: let I be any straight line with origin inO £ K™ 
and passing through the positive orthant then the point I D dR(X) is componentwise greater 
then I n dR(Y). If in particular both R(X) and R(Y) admit an ideal point x,y £ K™ then x 
is componentwise greater than y. This property describes exactly the idea that X is riskier 
than Y in terms of the geometry of the frontier of the sets R(X) and R(Y). A similar 
observation might be repeated for all the other properties: convexity/ sublinear ity state that 
the diversification/ aggregation decrease the risk of the portfolio and positively homogeneity 
excludes any liquidity risk. 

2. It is clear that if R(X) + K n = R(X) for every X £ Lf the region which lies above 
the frontier dR(X) represents higher level of risk that might be taken into account. This 
become a key point if we understand each component of R{X) as different criteria. Consider 
for instance two agents endowed with two different risk measures p\ and p 2 : we define 
R{X) = C(pi(X), p 2 (X)) where C{p\{X) 1 p 2 {X)) is a pointed cone in IR 2 . We will be able to 
find an agreement between the agents only if C(pi(X), p 2 (X)) n {(x,y) £ M 2 | x = y} ^ 0. 

3. Cash additivity in this multi- criteria framework slightly differs from the definition adopted 
in \22f . The interpretation is straightforward: if we add to the j th component of the portfolio 
X £ Lf a sure amount of money x then all the different criteria will agree that the risk 
decreases by —x. Every cash additive (d,n)-risk measure can be characterized by means of 
an acceptance set, namely a closed convex cone A C Lf and containing Lf(K). In fact the 
set valued map defined as 

R A (X) = l(Y^Xj)l n £ R n | x = (xi,...,x d ) £ M. d andX + x £ a\ 

is cash additive (d,n)-risk measure. Viceversa the acceptance set induced by a cash additive 
R is given by 

A R = {X G Lf | R(X) 4 n R(0)} 

Example 4 Definition^ show a simpler expression as soon as for every X the set R(X) satisfies 
R(X) + K n = R(X), as discussed in J22f and W^/ . On the other hand in certain cases we may 
need to provide a confidence interval of risk so that we expect the risk measure to take values in 
compact sets. We illustrate this idea through an easy example: again let 1, 2 be two agents endowed 
with two different risk measures p± and p 2 which are estimating the risk of a position X . Suppose 
that 1 only knows that 2 is more conservative (i.e. pi(X) < p 2 (X)), but does not know a priori 
the procedure used by 2. Thus agent 1 will provide a confidence interval [pi(X),M], where M 
represents the maximal capital requirement that 1 is willing to hold in order to cover the risk of 
X. On the other hand agent 2 will provide an interval [m, p 2 (X)] where m is the minimal amount 
of money 2 wants to save. Thus R(X) = [pi(X),M] x [m,p 2 (X)] describes the aggregate model: 
the agreement will take place only if R{X) fl {(x, y) £ M 2 | x = y} ^ which is equivalent to 
[ Pl (X),M]n[m,p 2 (X)]^0. 
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In the following Proposition we motivate the use of set valued risk measures instead of vector 
valued: if we consider different agents we are allowed to give different weights to each one of them, 
depending on the reliability. 

Proposition 5 let I be any straight line with origin in € R" and passing through the positive 
orthant R" . Consider the map p : U*/ — > R n defined as 

Pl (X)=ini{R(X)nl} (2) 

where the inf is to be intended componentwise and pi{X) = +ool„ if R{X) fl I = 0. If R 
is respectively monotone/convex/cash invariant/ 'sublinear /positive homogeneous then p is mono- 
tone/convex/cash invariant/ sublinear /positive homogeneous. 

Finally we state a well known automatic continuity result which is a key point for the opti- 
mization problems that we will consider in the core of this paper. 

Proposition 6 [22] Every (d,n)- coherent risk measure R : Uf =4 R d such that R{X) = R(X) + 
K n for every X £ , is continuous on . 



3 Optimal Portfolio Diversification 

Suppose that an agent is considering a vector X = (X\, ...,Xd) £ of risky financial positions. 
In a one-period model this vector might be composed by a basket of bonds, stocks, options so that 
every Xj > represents the value of the the i th position at the final time. This is the case of the 
two examples given in Section [5.1| For sake of simplicity we assume that the price of each asset at 
time is equal to ir(Xi) = 1 for i = l,...,d. Here the pricing rule ir is given endogenously, in the 
sense that prices are fixed a priori by market itself. The initial endowment x £ R will be given 
by a particular combination a £ R^ so that 1 = dtiir(Xi) = ^ i on. 

A different point (as the one followed in Section [5.2| is to consider a vector (Xi, X c i) that 
describes the possible losses/gains that the decision maker (as an insurance company) has to face 
holding the position JQ. 

In both cases an admissible risk diversification strategy will be given by any vector a £ R+ such 
that a i = 1) which represents the proportion of capital invested in each risky position. 
The decision maker is interested in redistribute and minimize his risk, by means of an optimal 
strategy. Namely given a real valued risk measure (as introduced in [10]) p : L°° — > R and X £ R d 
we have the following optimization problem 



min p (a ■ X) where A 1 

aGA d 




Here we extend the optimization problem ([3]) to the case of set-valued risk measures giving a 
different interpretation than the one in |22j , as explained in the following. Given a set- valued risk 
measures R : Lf =t R™ and X £ L^°, the agent deals with the following set-valued optimization 
problem 



vain R(aiXi, ...,adXd) where A' 

o6A d 




Since the vector X £ R d is supposed to be fixed we will consider, with a slight abuse of notation, 

R : A d =4 R", so that we will often write R(a) instead of R (aiXi, ctdXd). 

We suppose R" being ordered by the usual Pareto cone R™ which means a > b if and only if 

a — b £ R™ for all a,b £ R™. A pair (d,y), with y £ R(a) is an optimal solution to Q if 

R(a) C y - (R"\{0}) c for all a £ A d (see ;9» for more details). 

In the sequel we will consider the case in which the set R(a) takes the form 

R (aiXi, ...,a k X k ) := R(aiXi, a k X k ) + D(ai,a 2 , ■ ■ ■ , a n ), 
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Figure 1: A set- valued risk measure 



where € D(a) C R™ for all a e R+. 

A particular case happens when D = which leads to i? (a) = R (a) + . For instance when 
n = 1, the above condition degenerates to i? (a) = [R (a) , +oo). 

The vector R can be thus seen as a conglomerate of these different risk attitudes: some may be 
coming from regulatory agencies (external risk measures), others may describe different points of 
view arising inside the institution (internal risk measures). Risk managers would like to choose a 
strategy that minimizes all these different points of view, but this clearly becomes a quite tough 
task. 

Whenever the risk measure R : Uf =4 R™ is convex then R : A d =4 M." inherits convexity for any 
fixed basket X € R d . Thus the above program Q admits a solution since A d is a compact set. 
Moreover Q can be solved by searching for solutions to the following problem 

min R (o-iXi, ...,a d X d ) where A d =: \ a € Mi. \ V a, = 1 I . (5) 
« eAd I Y J 

In fact, it is easy to prove that solutions to ([5| are actually solutions to Q. Suppose that a 
is an optimal solution to ([5|. Since d is an optimal solution to ^ then it holds 

R(a) C i?(d) - (M;\{0}) c (6) 

for all a € A d . By easy computations we get 

R{a) = R(a) + D(a) C i?(d) - (M^\{0}) c + D(a) = R{a) - (M^\{0}) c (7) 

which shows that the pair (a,R(a)) solves Q. In the next section we propose a GP model for 
finding approximate solutions to (j5|. 



4 Risk Management through a GP model with satisfaction 
function 

In classical multi-criteria decision aid (MCDA) the agent has to consider several conflicting and 
incommensurable objectives or attributes which have to be optimized simultaneously. If D is the 
set of feasible solutions and fi represents the i-th objective function then the general formulation 
of a MCDA model is as follows [55]: Maximize (fi{x),f2(x),...f n (x)) subject to the condition 
that x € D. We suppose that each fa is continuous and D is a compact set which guarantee 
that Weierstrass theorem applies providing the existence of a solution. The Goal Programming 
model is a well known strategy for solving MCDA models; in this context, the agent seeks the 



G 



best compromise between the achievement levels fi{x) and the aspiration levels or goals gi by 
minimizing the absolute deviations (see [21 [201 [Ml HS1 HZ] ) and it can formulated as follows: 



in £ + (8) 



mm 

i=l 



t f + =g h i = l...n (Q) 

\ x e D y ' 

The GP model with weights (WGP), which represents an extension of ([8]), reads as: 

n 

min W i 6 t + W i S i ( 10 ) 
i=l 

s.t. / m + *t-^=9u i = l...n 

\ x e D y ' 

where wf and wf are weights or scaling factors. The GP model and its extensions have obtained 
a lot of popularity and attention because they represent simple models to be analyzed and imple- 
mented. However it is worth underlying that an optimal solution to {[sj) and ( 10 1 is an optimal 



solution to ([5]) if some optimality tests are satisfied ([M])- Among several extensions of these 
models which are currently available, it is worth mentioning the one developed by Martel and 
Aouni |26j which explicitly incorporates the agent's preferences. They introduced the concept of 
satisfaction function in the GP model where the agent can explicitly express his/her preferences for 
any deviation between the achievement and the aspiration level of each objective. In general, given 
three positive numbers £j , £<j and which will be called, respectively, the indifference threshold, 
the dissatisfaction threshold and the veto threshold in the sequel, a satisfaction function is a map 
F : [(),£„] — > [0, 1] which satisfies the following properties: 

• F(x) = 1, for all x € [0, where is the indifference threshold, 

• F(x) = for all x > where ^ is the dissatisfaction threshold, 

• F is continuous and descreasing. 

Depending on the thresholds' values, which strictly depend on the agent preferences, it could 
happen that positive and the negative deviations are penalized in a different manner. The GP 
model with satisfaction function is formulated as follows: 

n 

min^ W +F (5+) + wfF(Sf) (12) 
i=i 

fi(x) + 6f - dr = gi , i = l...n 
s.t. { x G D (13) 
5f,Sr g [0,&] 



Let us notice that the GP model ( 12 1 admits a solution because of the continuity of /, and F and 



the compactness of D. Some applications of this model can also be found in [H [H O [71 [S] . 

Let us now formulate a GP model with satisfaction function for risk management and optimal 
portfolio diversification based on the above multi-criteria optimization model ([5]): 



max]T w+F(Sf) + wfF(5f) (14) 



i=l 



Ri (a) + 5f - 5 i = g h i = l...n 
Ei on = 1 
oti > 0, i = 1 . . 
S+,Sfe[0,Cv} 



^ < %>o7i=i... P (15) 
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where gi represents the target of Ri for all i = 1 . . . n. If a solves ( 14 1 and a is an optimal 
solution to ([5]) then (a, R(a)) solves For this purpose one can use some of the optimality tests 
presented in [24]. 



5 Examples 

In this section we provide three possible applications of set valued risk measures to optimal risk 
diversification in presence of ambiguity concerning the risk criterion that should be adopted. In 
order to guarantee a simple computational procedure we choose a common mathematical frame- 
work for these examples, which is described in the following paragraph. 

In the first and second example we suppose that the agent is uncertain on the probabilistic model he 
will have to choose: namely a set of possible probabilistic scenarios {Qi}™^ 1 ^ s taken into account. 
Notice that Q x is excluded since it corresponds to the real probability P. Moreover we assume 
Qi << P so that even though P is unknown we have that the null sets are fixed a priori. The 
standard approach would be to define the coherent risk measure p(X) — sup i=2 „ +1 Eq i [—X] 
and solve the optimization problem given by equation (|3| . Here we will compare this deeply con- 
servative approach with a multi objective goal programming method. 

In the third example the reference probability is assumed to be known and we consider an agent 
that adopts the V@R as a risk estimator. Since we are considering a problem of diversification 
over risky positions and the V@R fails to be convex, the agent is uncertain on the criterium he will 
have to choose between the convex combination of the different Value At Risks J^i o-iV@R\{Xi) 
or the Value At Risk of the convex combination V@R\(J2i ctiXi). 

5.1 Robust methods for risk evaluation under model uncertainty. 
Illustrative setting In all the following computational examples we fix 

(n,^,P) = ([0 J l] 1 B [0 ,i],Le6) 

and the sequence of functions 

= -71 Vt = 2,...,n+ 1. 

Jo 

This last equation defines a sequence of probability Qi such that = fi and Q, ~ P. We simply 
compute J Q i^^duo = ^ (l — i) so that we deduce that lim^ /,(0) = +00 and lim^ fi(u>) = for 
every ui > 0. 

The linear case We firstly consider an example in which the agent has different beliefs in terms 
of probabilistic models but no risk aversion. We consider a general portfolio (a%Xi, adXd) € 
where a € A d . We define the set valued map R : If? =t R™ as 

n+l 

R{aiXi, ...,a d X d ) = [£%[- a ■ Z],+oo), 

where a ■ X is the usual scalar product in WL d . Simple computations show that R is a id, n)-risk 
measure satisfying (i), (ii), (iv), (v) and (vi) in Definition [2] 

In our illustrative setting we fix a portfolio composed by a non risky asset Xx(ui) = 1 and two risky 
assets X2(uj) = 2cj, ^3(0;) = 3cj 2 . The portfolio selection will be thus given by ai , aa , 0:3 £ [0, 1] 
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such that a\ + a 2 + a 3 = 1. We observe that 



Eq i [a ■ X] = Eq t [ceiXi + a 2 X 2 + 0.3X3} 



„ t-^duj Jo 



t (ai + 2a 2 oj + 3a 3 uj 2 )i u duj 

Jo 

2a 2 1 2a 2 3a 3 



6a 3 1 3a 3 



/ 1 1 \ / 1 1 2 

= ai + 2a 2 — - - — - + 3a 3 — - - — - - - — — — 
\lni i — lj \mi i — l (1— l)mi 

Notice that the classical approach would suggest to compute 

p n (a-X)= sup Eq i [—a-X] = — inf Eq^a-X] 

i=2,...,n+l t=2,...,n+l 

Clearly [a • X] > and Eq i [a ■ X] | a\ as i — > +00 so that 

p n (a • X) -ai. 

and the only strategy allowed becomes ot\ = 1. This means that the higher is the number 
of probabilistic scenarios taken into account, the more the decision maker will concentrate his 
investments on the non risky asset. 

We restrict ourselves to the case i = 2,3, 4. If we minimize Eq i [—a ■ X] separately we obtain 
the three possible ideal goals g\ = g 2 = 53 = — 1 with respect to the three different criteria. These 
goals are always realized by the strategy ot\ = 1, a 2 = a 3 = 0. In this way no risk is hazarded 
and consequently no real gains are reached. On the other hand if we choose three goals gi,g 2 ,g 3 
which allow a slightly higher level of risk then a better performing strategy arises. 
Let us choose the satisfaction function F(S) = ^g 2 , £„ = 60, and = u;j~ = = w 2 = 
I. We fix the goals g\, g 2 , and 33 to be equal to —0.7,-0.6,-0.5. The above GP 



U>q = III, 
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model for optimal risk diversification can be formulated in this setting as follows: 



max^K+F((5+)+ W rF(«5r)) 



(16) 

i=2 

(17) 

Eq 2 [-otiXi - a 2 X 2 - a 3 X 3 ] + 5% - 5 2 = g 2 , 
E Q 3 h a i x i - a 2X 2 - a 3 X 3 ] +6£ - 63 = g 3 , 
s.t < E Qi [-a 1 X 1 - a 2 X 2 - a 3 X 3 ] + 5+ - $7 = g iy (18) 
ol\ + a 2 + a 3 < 1, a\, a 2 , a 3 > 
„ < 5^,62,6^ ,5 3 ,6^,64 < £„ 

The numerical solution is performed by using LINGO 12 and provides the following solutions 
ai = 0.000000, a 2 = 0.9510359, and a 3 = 0.4896409S - 01. The decision maker clearly prefers to 
invest in X 2 even though the probability Qi(X 2 > Xi) < | for every i = 2, 3, 4. 
With respect to the classical optimal portfolio diversification in which the agent is forced to invest 
only on the non- risky asset, here we obtain a less conservative optimal solution which minimizes 
the distance between each achievement level and its goal. 

Entropic risk measure under ambiguity of the risk aversion. Again we suppose that the 
decision maker is uncertain on the probabilistic model Qj << P and his preferences are described 
by some exponential utility u(x) = 1 — e~ Xx . Moreover the agent might be more confident about 
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some probabilistic scenarios, so that his risk aversion A > will depend on Qj (i.e. A(Q,) = Ai). 
As usual the exponential utility induces the entropic risk measures 

p Qi (X) = mi{m e K | E Qi [u{X + m)] > u(0)} = — ln£ Qi [e~ A * x ] . 

A, 

We consider a portfolio {a\X 1: a^X^) S where a € A d and define a map R : =t W 1 as 

fe 

i?(aiXi, ...,a d X d ) = [pQ 4 (a • X), +00) 

i=l 

Simple computations show that i? is a (d,n) convex cash additive risk measure (i.e. satisfies (i), 
(ii), (iii) and (iv) in Definition [2]) . 

In our illustrative setting we fix a portfolio composed by a non risky asset X\(u)) = 1 and two 
risky assets X 2 (w) = 2uj, X 3 (uj) = (4w — l)l^ u> iy. The portfolio selection will be thus given by 
ai, 0C2, a 3 € [0, 1] such that a\ + a 2 + a 3 — 1. We thus have 

1 f Jq exp{-Aj(ai + 2a 2 w + a 3 (4w - l)lr u >ii) - ubii}du\ 
PQ z (a-X) = — In 1 1 



Computing explicitly the integrals one gets 

) oq i (a • X) = i In (i Ini) — i In (i — 1) — A^ai + 
Aj A, 



In 



' I _ e -|(2ct2Aj+lni) e -i(2c t2 A i +4Q 3 Ai+lni) _ e -(2a 3 A i +4a 3 A i +ln i) ' 

2a 2 A 4 + lni ' e" 3 (2a 2 Ai + 4a 3 A l + lni) 



Again we solve the optimal portfolio diversification problem via a GP model. We restrict ourselves 
to the case i = 2, 3, 4. If we minimize pq ; (a • X) separately we obtain the three possible ideal goals 
gx = -4.405194, g 2 = -5.497588, g 3 = -6.575445 with respect to the three different criteria. 
As in the previous model, let us choose the satisfaction function F(S) — 001*52 , £ v — 60, and the 
weights wf — = = w 2 = w t = w 3 = 5- We set the parameters A 2 , A3, and A4 equal to 4, 
5 and 6, and the goals g±, g 2 , and 173 equal to —4,-5, and —6 respectively. The above GP model 
for optimal risk diversification can be formulated in this setting as follows: 



s.t. 



(w+F(S+) + wrF(S-)) 



Pq 2 (aiXi + a 2 + X 2 + a 3 X 3 ) + 5 2 - 5 2 
p Qa (aiXi + a 2 + X 2 + CK3X3) + 8$ - 83 
PQi {ol\Xi + a 2 + X 2 + a 3 X 3 ) +64-64 
ai + a 2 + a 3 < 1, at\, a 2 , a 3 > 
{ 0<6+, 6^,6+, 6^,6+, 64 <& 



92, 
93, 
34, 



(19) 
(20) 

(21) 



The numerical solution is performed by using LINGO 12 and provides the following solutions 
ai = 0.6659041, a 2 = 0.000000, and a 3 = 0.3340959. 



5.2 A new point of view concerning V@R 

In the previous examples we built up a set valued risk measure aggregating different real valued 
risk measures that were generated different probability beliefs. Moreover in both example R 
depended on the vector X = (Xi, X^) 6 U£ only through the sum of the components. In this 
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last example we suppose that the historical probability measure is known: in this framework the 
most popular (and also most debated) risk measure is V@R\ defined as 



V@R X (Y) = - sup{m e K | ¥(Y < m) < A}, 

where A € [0,01,0.05] and Y € L° is any J 7 - measurable random variable. Notice that for every 
a > we have V@R\{aY) = aV@R x (Y) and V@R\(Y + c) = V@R X (Y) - c for every c € R. 
Nevertheless the Value at Risk is not convex on the space of random variables and for this reason it 
does not sense the effect of diversification. This lack has an immediate consequence: if we consider 
a basket of financial instruments X = (Xi, ...,X d ) e Uf and a e A d we cannot guarantee any 
order relation between 



V x (a) ~J2^V@Rx(Xi) and V x (a) =: V@R X PTa^ (22) 




As a consequence the decision maker should be uncertain among all the possible values x € 
[V\(a) A V x {a) 1 V\(a) V F A (a)]. This problem can be clearly reinterpreted via a multi objective 
goal programming. 

Another common risk measure is given by p w (Y) = — ess inf F with Y € L°°, which is a coherent 
(and thus convex) risk measure. The clear drawback is that worst case risk measure is too restric- 
tive and conservative from the point of view of an agent who is investing his capital. Anyway it 
can be exploited to give an upper boundary of the maximal capital requirement necessary to cover 
any possible expected loss. 

We introduce the following set valued map R : R 2 defined as 



R(X u ...,X d )= { 



V x (l d ) < x < Pw (l d ) 
(x, y) e R 2 : and 

v\u) < y < vHu) + Z{\fXll d J (» - VxiU)) J 



(23) 



where V\(l d ) =: E 4 V@R x (Xi), V x (l d ) =: V@R\ (J2i X t) and p w (l d ) = E t P W (X t ). Notice that 
the set R(Xi, X d ) is a triangle and has a minimizer and a maximizer (w.r.t the Pareto cone 
K+) given respectively by (V\, V x ) and (p w ,p w ). 

We consider a vector (aiXi, ...,a d X d ) e where a G A d so that by the positive homogeneity 
of V@R and p w we find 



R{a\X\, a d X d ) = < 



V\(a) <x< p w (a) 
(x, y) e K 2 : and 

V\a) < y < V x {a) + ^{tXtt ( x ^(«)) 



where p w (a) = J2i a iPw(Xi). 

As usual our illustrative setting is given by (f^J 7 , P) = ([0, 1], #[o,i] , Leb) but in this case the 
Lebesgue measure in chosen as the reference probability. In order to clarify the example we 
consider three financial positions which allow negative losses namely X\{ui) = 0, X^uf) = 2u> — 1 
and X 3 (uj) = 2 — 4w 2 . Since the three random variables are continuous we deduce that for A = 0.05 

V@R\{X 1 ) = V@Rx(X 2 ) = 0, 9 V@R X (X 3 ) = 1, 61 
PviXi) = Pw (X 2 ) = 1 Pw (X 3 ) = 2 

so that Vx(ct) — a>2 ■ 0, 9 + a 3 ■ 1, 61 and p w {a) — a 2 ■ 1 + a 3 ■ 2. We need to compute V x (a): notice 
that Y,i a i x i(u) = a 2 (2w-l) + a 3 (2-4w 2 ). In particular £\ ^^(0) = 2a 3 -a 2 = -^ i a i X i (l). 
Then the function J^^ajX^a;) have a maximum point in ui = f^. In general the minimum of 
the parabola /(w) = a 2 (2uj — 1) + 03(2 — 4w 2 ) will fall on uj = if 2a 3 — a 2 > and on u> = 1 if 
2a 3 — a 2 < 0. 
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(1) Suppose that a 2 — 2a 3 . In this case the V@R can be simply computed as 

V@R X (a • A) = f Q 

(2) Assume a 2 < 2a 3 . We find that the solution of /(0) = f(uj) is given by uj — 0; f 3 ^- and 
P (a ■ X < /(0)) = 1 - We have two possible cases 

(a) A < 1 - |f- then 

F@i? A (a ■ X) = f (1 - A) . 

(b) A > 1 - %l th en 

V@R X (a-X) = f 



(3) Assume a 2 > 2a 3 . We find that the solution of /(l) = f(uo) is given by co = 1; ^ — 1 and 
P (a • A < /(l)) = ^ — 1. We have two possible cases 

(a) A < g - 1 then 

y@i? A (a • A) = / (A) . 

(b) A > ^ - 1 then 

(\ _ (_2^_ — l)\ 
1 2 J ' 

Finally, fixing A = 0.05 we can conclude 

if a 2 < l-9a 3 then V@R X {a ■ X) = 0.9a 2 - 1, 61a 3 

/ \ 2 2 

if 1.9a 3 <a 2 < 2a 3 then V@R X (a ■ A) = -4a 3 ( - 0.475 ) + p- + 2a 3 - 1.95a 2 

\4a 3 J 2a 3 

if 2a 3 <a 2 < 2.1a 3 then V@i? A (a ■ X) = -Q.9a 2 + 2.99a 3 



if2.1a 3 <a 2 then V@R\ (a ■ A) = -4a 3 — ^ + 0.475 + + 2a 3 - 0.05a 2 

V 4a 3 / 2a 3 

The above GP model for optimal risk sharing can be formulated in this setting as follows: 

2 

max]T (w+F(S+) + wr F (Sr)) (24) 

i=l 

(25) 



s.t. < 



f E-=i «i^@i? A (Aj) + St - Si = g u 
V@Rx(J2ti <XiXi) + &t - S 2 = 92, (26) 

cti + a 2 + a 3 < 1, ai,a 2 ,a 3 > 
{ 0<5t,S^,St,d 2 <£ v 



We choose the goals g\ and g 2 to be equal to 0.3 and g 2 = 0.5 respectively. LINGO 12 provides 
the following optimal solution: a x = 0.6823770, a 2 = 0.2151639, and a 3 = 0.1024590. 

6 Conclusions and further developments 

The recent notion of set-valued risk measure appears as powerful tool that can be exploited 
to overcome many complications that arise in risk management. Risk, understood as capital 
requirements needed to cover expected future losses, becomes an ambiguous factor to determine 
as far as a manager has to face different criteria or is uncertain on the real probabilistic model that 
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lays underneath the financial problem. Inspired from the statistical notion of confidence intervals, 
the general definition presented in this paper allows to consider compact valued risk measures. In 
this way we associate to any financial position a cloud of different risk levels, instead of a single 
number, taking into account all the multiplicity of ingredients that characterize this computation. 
As illustrated by some examples we may thus formulate an optimal risk diversification problem 
which allocates the risk of a given portfolio in an optimal manner. This is a set valued program 
which can be reduced to a vector valued model if the images of the set valued mapping admit an 
ideal point. Using a Goal Programming approach with satisfaction function we are able to provide 
approximate solutions to this vector model: the presence of several free parameters is the strength 
of this approach since this allows a calibration of the model sensitive to the risk aversion of the 
agent. 

We have then illustrated three different examples which support this approach: in the first and 
the second one, the agent has fixed a risk procedure p but he is uncertain about the probabilistic 
model Q. In such a case the functional form of pq(-) will explicitly depend on Q and the standard 
literature would suggest to take a supremum suPqPq(-) to compute the capital requirement. As 
explained above, such a strategy would often force the agent to avoid any risk in his decision. 
Through the GP model we find a non trivial diversification strategy which takes into account all 
these different possible scenarios Q. In the third example we provide a case of compact-valued risk 
measure built up from the celebrated Value at Risk. As well known, the Valued at Risk is convex 
only on the space of Gaussian random variables, but it looses this property if we extend its domain 
to more general random variables. As a consequence the V@R is not sensitive to diversification 
and for this reason it might not fit our optimization problem. The method we have proposed in an 
illustrative setting is a natural starting point to overcome this controversial and debated feature 
of the V@R. 

For future developments, we are going to conduct a statistical analysis of the model illustrated 
in this manuscript by using real data and by estimating the images of a set-valued risk measure 
through the analysis of confidence intervals. 
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